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ABSTRACT 

This  report  presents  the  theoretical  and  experimental  results  related  to  flows  over 
such  simple  bodies  as  semishpere,  cone,  and  wedge.  These  were  investigated  in 
hypersonic  wind  tunnels  of  various  classes  at  nearly  equal  Mach  and  Reynolds 
numbers  (M  from  7.0  to  8.0,  Re0  from  135  to  240)  but  at  notably  differing 
free-stream  velocities:  790,  2700  and  6000  m/s. 

Hypersonic  resistance-  and  arc-heating  wind  tunnels  and  a  hypersonic 
MHD-accelerator  wind  tunnel,  as  well  as  their  test  equipment  are  described.  A 
technique  to  determine  parameters  of  a  free  stream  over  a  model  including 
gasdynamic  ones  (M,  Re,  Re0,  T,  P,  V,  Ho,  H)  and  physico-chemical  ones  (y,  T0,  Tv, 
Tn,  Q),  as  well  as  a  verification  procedure  are  presented. 

Experimental  data  on  pressure  distribution  over  models,  shock  wave 
positions/shapes  were  obtained.  The  experimental  data  are  compared  with 
respective  calculated  results  obtained  by  using  the  VSL  theory  and  the 
Navier-Stokes  equations  for  the  test  conditions  adopted.  Interrelation  of  the 
theoretical  and  test  data  is  shown  not  to  be  unambiguous,  especially  in  what 
concerns  locations  and  shapes  of  shock  waves. 

The  VSL  theory  is  employed  to  compute  local  parameters  of  the  shock  layer.  The 
influence  of  alkali  metal  seed  on  the  rate  of  relaxation  in  the  shock-layer  is 
considered.  Features  of  profiles  of  both  temperature  and  air  component 
concentrations  in  the  shock  layer  are  shown  to  be  useful  in  explaining  a  set  of 
phenomena  revealed  in  experiments:  a  radiation  offset  upstream  of  a  sphere 
hypervelocity  shock  layer,  a  radiation  relaxation  at  a  cone  apex.  The  impact  of  the 
seed  metal  is  shown  to  be  insignificant. 


INTRODUCTION 


The  development  of  a  new  generation  of  hypersonic  transatmosphenc  vehicles 
featuring  high  L/D-ratios  such  as  NASP  requires  that  the  physico-chemical 
processes  inherent  in  flows  over  vehicles  be  studied  more  comprehensively,  on  the 
one  hand,  and  such  test  facilities  be  constructed  which  are  capable  of  reproducing 
these  processes  adequately,  on  the  other  hand. 

The  ranqes  of  speeds  and  pressures  for  which  these  studies  should  be  earned  on 
are  governed  by  flight  paths  of  such  vehicles.  Their  acceleration  to  orbital 
velocities  occurs  in  relatively  dense  atmosphere;  they  can  also  perform  maneuvers 
at  altitudes  of  50  km  to  70  km  and  at  flight  speeds  of  M  =  20  to  M-25. 

Many  known  articles  are  devoted  to  the  description  of  flow  fields  of  chemically 
reacting  gases  both  over  simple  bodies  and  real  vehicles  [1],  |2J,  PJ- 
The  thermochemical  models  underlying  these  calculations  and  the  numerical 
algorithms  developed  by  different  authors  differ  considerably  because  of  a  number 
of  unverified  assumptions. 

In  the  calculations,  the  differences  in  reaction  rate  constants  for  flows  in 
thermodynamic  equilibrium  are  not  as  important  as  nonequilibnum  ones  ln 
non-equilibrium,  the  flow  field  depends  greatly  on  the  thermochemical  model  used 
Sy  the  authors.  In  most  previous  investigations,  the  reaction  rate  is  supposed  to 
depend  on  translational  temperature  only.  In  latest  publications  [4],  [  1- 
chemical  reaction  rate  is  assumed  to  depend  both  on  translational  and  vibrationa 
tpmneratures  because  of  a  close  relationship  of  dissociation  with  vibrational 
dearees  of  freedom.  The  temperatures  of  rotational  and  electron-vibrational 
decrees  of  freedom  are  also  assumed  to  be  different.  Considerable  gradients  in 
concentrations  of  different  air  components  and  excited  electron  states  occur  under 
the  flow  conditions  behind  the  shock  wave  front.  The  tests  earned  out  in  shock 
wind  tunnels  [61  show  the  concentration  of  NO  to  vary  most  rapidly  behind 
Zck  wavX  ^Jh  this  concentration  starting  to  rise  immediatelly  following  the 
dissociation  A  rapid  NO  electronic  excitation  behind  the  shock  wave, 
ST<o“h  transnational  temperature,  leads  to  an  Intensive  ul.r.vtolet 

emission  in  the  Gamma  band  system  (>.-200  to  320  nm). 

The  situation  is  further  complicated  by  nonequilibrium  emission  caused  also  by 
otlTe^air  components  both  hi  visual  and  ultraviolet  and  infrared  regions.  The 
intensity  of  nonequilibrium  emission  can  exceed  that  of  equilibrium  em^°n  ^ 
factor  of  2  to  15  [7].  However,  these  calculations  [8]  are  not  confirmed  by  flight 
Jest  results,  although  the  intensity  of  radiation  heating  exceeds  convective  heating 

intensity. 

All  these  considerations  dictate  that  the  thermochemical  models  and  mathematical 
XorithL  be  trifled  in  laboratory  conditions.  It  is  desirable  to  choose  those 
sneed  and  pressure  conditions  which  could  allow  comparison  of  the  results 
obtained  by  different  authors  in  various  facilities  and  of  respective  numencal 

algorithms. 

The  simplest  technique  to  reproduce  the  conditions  of  hypersonic  equilibnum 

l.ln'Tcin.9^.*  to  reproduce 
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different  properties  of  the  test  gas  (F™on-U  h-  ^  ^ 

(y  =  1.67))  over  conical  bodies  are  mves  ig  ^  stand-off  distances  and  shock 
some  distributions  of  surface  pressures,  cone  angles 

shapes.  It  is  “erred  from  these :  data  “  real  g.s^  ^  ^9  ^  ^ 

ess 

Scfe  lVSniqSor  their  investigation  is  given  in  (111  which  also  contains  a 
description  of  their  merits  and  weaknesses. 

The  authors  believe  that  such ,  facilities 

shock  wind  tunnels  and  shoe  win  been  considered.  Among  these,  optical 

experiments  carried  out  in  th  determination  of  shock  stand-off  distances 

S^pes^afStfenges  oAheir  possible  application  depending  on  gas 

ZZ  consideration  is  given  to  £  — 

^“fedS'riSoS shock 

Another  type  of  test  .acuities 

V-6  to  8  km/s  is  ballistic  r“9e ffiq  /JoC  Research  Center  aimed  at 

SLOT*  -  vaSf.  models  ^tS.TcSp^"sh%S 

-  '“*u“s  - 0,6  COT,OT 

Sto  ^°lds  — S  redistributing  heat  flukes 
^eofeLTlnd  experimental  data  on 

sphere,  a  wedge,  a  corns  • j  ^d“  “d„ lf^es  .1  4  to  8  km/s  and  pressures 
PAET  are  compared  m  [15]  tot  the  9  obtained  in  shock  wind  tunnels  and 

of  2- 10  to  104  Pa.  All  the  experimental  data  b  ^  di{ferent  thermodynamic  gas 

ballistic  range  facilities.  The  comp  .  erature  model  with  seven  chemical  air 
models.  It  is  conclude  a  t  with  experimental  data  than  e 

components  is  in  better  agreemen  rfect  gas  and  a  thermodynamic 

one-temperature  model.  The  model^  P"^  dJcussion.  Unfortunately,  the 

equilibrium  gas  is  inadequa  Depending  on  the  level  of  model  surface 

test  model  material  used  is  no  P  m  vary  under  such  conditions  by  10 /o  to 
cataheity,  the  shock  stand-off  distance  may  va  y 

18%  [16]. 
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Short  run  times  in  shock  wind  tunnels  and  small  dimensions  of  models  in  ballistic 
test  rigs  restrict  the  capabilities  of  experimental  investigations  of  nonequilibnum 
processes  considerably. 

The  application  of  classical  wind  tunnels  for  these  purposes  is  impeded  because  of 
technical  and  physical  reasons. 

It  is  shown  in  [171  that  in  order  to  reproduce  hypersonic  flight  conditions  in  a  wind 
tunnel  test  section,  three  conditions  must  be  satisfied  simultaneously: 

—  >H;  ^->Edis+Ev; 

2  ~ 

xchem.m  ^  xchem,n_  . 

xaerod.m  xaerod.n 

As  for  classical  wind  tunnels,  they  prove  to  be  incapable  of  realizing  in  M 
the  first  or  the  second  condition.  Even  if  it  could  be  possible  to  overcome  technical 
problems  of  achieving  and  maintaining  the  plenum  chamber  gs  parameers 

required  to  obtain  a  test  section  flow  with  full-scale  velo^es^ep  u>  3  Eclated 
the  studv  of  nonequilibrium  processes  would  be  impossible.  When  the  dissociated 
nas  is  expanded  in  a  hypersonic  nozzle,  a  considerable  portion  of  stagnation 
^  i  becomes  "frozen”  in  internal  degrees  of  freedom  and  chemical  reactions. 

“  »  amount  to  50%  and  more.  The  specific  heat  ratio  is 

also  "frozen",  and  the  flow  over  a  body  becomes  equivalent  to  an  ideal  gas  flow. 

A  hvnersonic  MHD-air  flow  acceleration  wind  tunnel  constructed  and  operated  at 

ftots? 

Reynolds  numbers,  but  ydth  significantly  different  rffects  of  P^*1 

Kr/cel  r„e,  K  c^SiJ’Sg  - 

22%  f  pp-f  E— ; 

fo  be  rather  effective  in  the  research  of  gas  dynamics  and  heat  transfer  to  scram, et 
combustion  chamber  components. 

^“Sd'f^hSr^retSS,  results  obtained  by  other  investigators. 

already  stuaiea  y  t  rfil  ancj  rin  the  flow  velocities  required 

££  biUTSf  t„A“Ss9.nd  flol  pressures  horn  40J>. f»  «0 
parameters  are  also  characteristic  of  an  intensive  nonequilibnum  radiation 
shock  layer. 

facilities  of  different  types. 
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DESCRIPTION  OF  WIND  TUNNELS,  INSTRUMENTATION  AND  MODELS 

Ohmic  (resistor-type)  heating  wind  tunnels  have  found  a  wide  application  in  the 
investiaatlonsof  hypersonic  perfect-gas  flows.  The  arc-heating  facilities  are  capable 
of  producing  flows  of  both  vibrationally-excited  gases  and 

Tests  were  conducted  at  stagnation  temperatures  of  400  K  to  1000  K  and  251)0 
4000  K  and  flow  velocities  of  800  m/s  to  1200  m/s  and  2500  m/s  to  3000  m/s, 

respectively. 

A  hypervelocity  MHD-gas  acceleration  wind  tunnel  makes  it  possible  to  achieve 
flow  velocities  of  4-103  to  8-l0»  m/s.  Since  the  values  viscosity  of  g0  ^e  ^her  great 
because  of  a  high  stagnation  enthalpy  past  a  shock  wave  near  a  tes*  , 

Revnolds  numbers,  Re0,  will  be  relatively  moderate.  Therefore,  the  Reynolds 
numbers,  Re0,  realized  in  this  wind  tunnel  is  considered  to  be  a  basis  from  which 
to  choose  other  wind  tunnels  and  their  operation  regimes. 

rs.  “ 

seems  to  be  impossible. 

In  view  of  the  above  considerations,  the  hypersonic  low-density  wind  tunnrf 
VAT-102  has  been  chosen  to  reproduce  this  velocity  range  l2^  ^XShig 
Mach  number  M=8  realized  in  this  wind  tunnel  serves  as  a  cntenon  for  ch  g 
Mach  numbers  required  for  the  present  investigations. 

both  nonequilibrium  and  equilibrium  flows  with  high-stability  parameters. 
Comparison  of  wind  tunnel  characteristics  in  terms  of  the  level  of  their  testsection 

bi°49^ 

t»db<o„  te  pampas.  aW  wave  <o,  ,ha  wiad  .uane.  VAT-.02  *1.  be 
characteristics  are  given  below. 

nozzle  raaiataine^  by 

instafled  at  their  exit.  . 

The  given  tests  were  conducted  using  a  contoured  nozzle  with  the  throat  size  of 
7  5  ?o-3  m  and  the  exit  diameter  of  0.185  m,  the  flow  core  diameter  being 
0  065  m.  Se  measurements  included  plenum  chamber  gas  pressures  and 
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fpmneratures  as  well  as  flow  fields  past  shock  waves  P^,  at  nozzle  exit  and  along 
iStraults  obtained  are  given  in  Figs.  2  and  3  .According  to  these 
measurements,  the  Mach  number  M=7.63  at  P„=2  Pa  and  respective mode 
_  _  u  -Q  —070  tn  'lOO  Re  «3800  to  4000  were  achieved  m  the  test 

ones.  The  measurement  error 

was  ±3%  for  pressures  and  ~1%  for  mean-mass  temperatures. 

Static  test  section  flow  temperatures  were  ~23  to  27  K  which  meant  that  they  were 
Static  test  section  u  p  rf  nitrogen  equilibrium  condensation  values.  But 

Ebd^r^Sn-“gtt^  "  P“  1 

condensation. 

ssS? Hr=£i=  SssSS j 

KSSSS'SSSs 

5*J£3SSH?SSS£ 

photographs  were  made. 
m  Omit  o'  Ite  tvind  tunnel  Wtlt  SiL 

tunnel  consists  oi  the  following  “^““Ibe”  J£h  widows  to  observe 

nozzle,  a  test  section  tn  the  form  of  an  E.ffel  chamber  wt  ^  ^  exhaust 

The  arc  heater  has  a  £p«slb.e 

to  m^m^^gh-ftegum^^dischargecunent  and^voltag^flucUl^ibM^nc^fi^^11^ 

~ 

of  electric  discharge  occurring  in  it  is  given  in  [23,24]. 

...  .  j  j  ikali  seeds  to  the  heater  flow  which  enables  vanous  optica 

m^g5mns',0„ddb^es.  Sees  it  possible  have  an  almost  equilibrium  flow. 

,  1  i  mm2  pYit  section  was  used  in  these  tests.  In 

a  number  o7oLT2sL,Tnozzle  based  on  the  secondary  MHD-accelerator  nozzle 

r.:— .  system  designed  >“  ^^^tTg^XSl^f 
and  alkali  seed  flow  rates,  pressures  heater ^current  ^  ■  as  well 

static  pressures  along  the  nozzle  and  static  w*sanamv*  exit. 

as  pressure  fields  past  »  Within  required  ranges 

Magneto-inductance  gages  wereuse^  t  ?unnei  was  mn  as  foUows.  After 

Sfind  mSrel  wa's  s?.Sli  id  the  required  test  regime  stabilized,  the  model  was 
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introduced  into  the  flow.  After  some  5  to  7  s  the  model  was  removed  from  the  test 
section  and  the  wind  tunnel  switched  off. 

in  some  tests,  .Umll  seeds  were  injected  into  the  flow;  this  was  done  when,  as  , 
rule,  the  model  was  already  in  the  flow. 

The  parameters  that  were  impossible  to  measure  at  the  nozzle  exit  were  determined 

using  known  g.s  dynamic  relations  lor  an  ^^/"g^fseSTon  positions, 
and  lor  a  flow  with  instirtijjnmly  dls, buttons, 

"  ”°“'e  ““ 

the  present  tests  are  given  in  Fig.  6.  .... 

”  £ wlnd  ,omel 

and  Fiqs.  8  and  9  its  general  views.  .  . 

The  operation  of  the  wind  tunnel  was  conducted  as  follows,  A  g.s  heated  nr  lm.te 
l tie  operation  ui  oin5p  -  qiinnlied  to  mixing  chamber  (Z) 

(1)  up  to  To~3500  K  at  0  o  a  o^  ht)  were  injected  by  metering  device  (3). 
where  slightly  ionized  seeds  (U  of  we >  g  )  Jersonic  *nozde  (4)  entered 

Then  the  gas  passing  throug  P  ^  constant  electric  and  magnetic 
MHD-accelerator  (5)  where  it  accelera  expanded  to  specified  static 

fields  up  to  required  velocities.  M  er  the  gas  was  p  section  (?)  ^ 

'^"h.re  through  a  dUhrser,  a  coohng 

device  and  a  set  of  ejectors. 

The  supply  of  seeds  in  the  form  of  KN.  eutectics  to  the  ^Hng  dScfS  w“ 
the  wind  tunnel  plenum  chamber  was  ^c^^aJ.ement  supply  system  whose  supply 
a  combination  of  a  plunger  press  and  optimal  placement  and  a  special  injector 
lines  were  protected  by  an  inert  ^  S"  distribution  across  flow 

ensured  finely-ground  seeds  (da  f  ™^  taJnducting  flow  was  formed  that  featured 

another  high  stability  of  its  parameters  [27]. 
SS^anTr^oSblts  outleT  parameters  were  matched  with  th. 

inlet  parameters  of  the  MHD-accelerator  channel.  chamber 

The  MHD-accelerator  consists  of  two _  mam  ^tS-  payday- type  rectangular 
(MHD-channel)  and  an  eJectrom  g  ^  walls  d  two  isolating  walls 

MHD-channel  includes  two  sectionahzed l  f^trohe  ^  to  B  =  2.5  T  is 

perpendicular  to  them.  The  magne  ic  ^  E -shaped  parts  with  an 

produced  by  an  electromagne  ^possible  to  apply  MHD-channels  of  different 
MHD-channel  between  its  bars.  “  P  .  Ltemdimensions  are  also  varied  as 

td^d^tiT oStS  design  features  is  9-^  ^ 
Those  tests  used  a  MHD-channel  550  ™  “  0’elSrede  pairs  and  elecuode 

produced  along  a  length  of  43 m  ^dfte  umtom  8^  within  the 

electrodes  were  used. 
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Tile  accelerator  power  supply  ™  ^ 

containing  1  to  20  isolated  secondary  tt “^hTSpi*.  bridge-circuit 
connected  to  the  respective  e  ec  P  h  ke  Thyristor  voltage  regulators 

S  d“  “  “,ive  io*ds  did 

rrJ—  witb  15.  10  and  ,0  power  circudtand  voH.ges  o,  000  V.  000  V 
and  900  V,  respectively,  were  used  in  the  presen  accelerator 

Removable  secondary  rectangular  nozzles  were  attached  to  «  wth  the 

channel  through  an  insulating  spacer.  The  ™ri**>*  ™ 

-  - — -  -  - 

Z  rs~  designed 

l“cTL"L°g ^“to  flow9  The  -«£—  -  te  “ 

,  in  the  hvuersonic  wind  tunnel  with  an  arc  gas  heat  . 

Te  wind  tunnel  SSTot^efwShTe 

35  -  S  V  nd  ^“'^"erraUeet^ 
simultaneous  photographing  an  P  seed  supply  was  stopped 

deenergised  and  the  modd ’ '^SeSlagnet  and  recording  devices  were 
+>»p  heater  was  switched  oil.  me 

Err rr^ssStf.tasss 

gas  and  seed  flow  rates,  pres®.  .  1,  ti  s  alonc,  MHD-channel  and  secondary 
electromagnet,  static  pr< sss™e< Sons  on  electrodes  along  MHD-channel  length, 
nozzle,  current  and  voltage  distabut  normal  shock  P6  in  test  section.  The 

static  pressure  and  pressure  fi  P  inductance  pressure  transducers  of 

differerfl  ^  ^  ^  ^  ^ 

During  the  preparation  f°r  u^d 'oscillograph^  was°Shed  over  to 

provide  for  a  more  rehabie  msula  f  P  lvanic  isolation  systems  m  the 

components  and  the  application  oi  cornea  ^/measurement  error  for  pressure 
circuits  used  to  xneasum  electric  was  no  more  than  2%. 

was  no  more  than  fa  f  MHD-cras  acceleration 

The  evaluation  of  flow  Para”®te^  “  ^®han  to  "conventional  wind  tunnels,  since 
facility  is  a  more  complicated  process  U 1  the  channei  under  the  action  of 

stagnation  parameters  vary  citation  Direct  measurements  of  gasdynamic 

pondermotor  forces  and  Joulefl/f  *P//°nt  considerable  difficulties.  To  facilitate 
parameters  in  a  high-energy  P  d  on  a  numerical  solution  of  the 

their  determination  a  Vr0C^ This  procedure  is  effective  when 
so-called  "inverse  flow  problem  applied  as  boundary  and  initial 

si  SKSS£^£^  - - — 
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The  MHD-channel  flow  parameters  not  subject  to  measurement  trore  calculated  in 

IZSSSEZJSSSTS  sssn.xss ?£? “  Z  — 

form: 

du  dP  j  o  or  Pu  • 
pu— =  -  — +  lyB  2Cf  r,  1 
dx  dx 

pu2  —  +  pu  =  jyEy  +  jxEx  -  Qw  '< 

H  dx  dx  * 

pUF  =  m;  H  =  H(P,T); 


T-fxl  s  Ve(x)-AVnear. 

lyW-^r-  h_.5T' 


•  -  r_1i  • 


r(n-H?i) 

r8PH 


Geometric  channel  parameters.  matured 

SeCSio^  F”rp“.mJ:Se™re9d  a.  accelerator  channel  input  were  used 
as  initial  conditions.  ...... 

oTthe  calcnUted  dam  ebtah,^  £  rite 

S"”  SioM  from  nm  to  mn  do  no. 

GX.m^pt«m.e.ers  and  Ur.rmodyn^c  toncttons  to  flre^ 

~$sssr£ -  p"“^  “  ob,aiMd 

Z"n,  ert-  -  £  « 

SLSSSSi  Z  test  ranges  of  temperatures  and  pressures: 

o2+m  =  o  +  o  +  m 
n2  +  m  =  n  +  n  +  m 
NO  +  M  =  N  +  +0  +  M 
N  +  02  =  NO  +  O 
0  +  N2  =  NO  +  N 
N  +  O  =  NO+  +  e 
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TOese  reactions  were  accomplished  by  ionisation  reactions  between  K  and  N. 
seeds: 

K  +  e  =  K++e  +  e 


4-  p  =  Na+  +  e  +  e 

Na  _  ,  T  KTO  n  N  NO+  K.  The  equations  were 

Here  M  is  any  particle  oi  O* ^TOrrSittal  oi  calculated  parameters  to  test 

U  L P^in"r  sLnda.  norsle  and  a.  its  -t 
accelerated*  irTtt^WHD-<m^nelthey  wme^'mclud^rir^er^e'^Se^  grant,  into 


equation  for  full  enthalpy  H  is: 

A(pH0uF)  =  ^kP. 

dx  D 


(i) 


re  coefficient  P  represents  complete  energy  losses  (including  energy  losses  for 
ation,  electrode  losses  and  heat  fluxes  to  walls). 


ation  electrode  losses  diiu 

equation  for  energy  flux  inherent  in  N2  molecule  vibrations  has  the 
fEVi  pjUpj  _  pt ~gUBh  ~Ifc  +  (ZlVT  +  ZlW'  +  Zlchem^  ’ F 


(2) 


“t'iKS  a»^  SSir^nTrr" 
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Rates  Ty  of  energy  transfer  to  vibrations  in  the  processes  VT,  W  and  in  chemical 
reactions  are  calculated  by  using  relations  given  in  [35  to  39]. 

The  equations  for  02  and  NO  molecule  vibrations  have  the  similar  form: 

—  f — — -  ^ l  =  ^2VT  "t  ^2W  "t  ^2chcm  ^ 


F  dx  l  m2 


J-—  5^-p3uF  ]  =  Z3vt  +  Z3VV'  +  Z3chem 


F  dx  ^  m3  ) 

The  gas  composition  is  determined,  accounting  for  the  following  reactions: 

N2  +M  <-»N  +  N  +  M 
02  +M  <-»0  +  O  +  M 

NO  +  M<-»N  +  O  +  M  (5) 

02  +  N  ONO  +  O 
N2  +  O  ONO  +  N 
N2+02O2N0 

The  equations  for  particle  flows  are  written  as  follows: 

iAfPLuF^/j  (i  =  1,6)  (6) 

F  dx  nij  ) 

Chemical  reaction  rate  constants  (6)  required  to  calculate  functions  of  particle 
collision  J,  are  taken  from  [35  to  37]. 

The  total  enthalpy  Ho  is  derived  from  relations  from  [35  to  37]: 

Hn  =^  +  H 


H=I£LHi 

i=l  P 

H;  =  RT-^  +  —  + 

1  m;  mi 

7t;  =  —  for  i  =  1,3; 


it,  =  —  for  i  =  4,6 
2 


The  equations  for  momentum  flow  are  used  in  the  form: 

A[(p  +  Pu’)p]-2C,££+if^pf  ™ 

The  set  of  equations  (1)  to  (10)  including  the  equations  for  mass  flow 

—  (puF)  =  0;  (9) 

dx 

and  for  state 

_VJR  T  (10) 

P  ^  mi  ‘ 

and  experimental  values  of  Vk(x),  Ik(x)  make  it  possible  to  calculate  gas  parameters 
in  the  channel. 
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The  calculation  of  the  parameter  p,  being  dependent  of  electric  field  strength  is 
impeded  by  the  field  nonuniformity  across  the  channel.  This  parameter  is  varied, 
herefore,  in  the  range  of  0  <  p,  <  1.  But  its  influence  proves  to  be  insignificaiit 
because  intensive  deactivation  of  N2  vibrations  during  V-T  relaxation  [40  to  43]  m 
the  case  of  collisions  with  N,  O  and  Na  atoms  results  in  equilibrium  vibrational 
and  translational  degrees  of  freedom  of  N2  molecules.  The  final  calculated  data  are 

given  at  p^O.5. 

The  parameters  p,«0.8  and  0*5-10-3  were  chosen  from  the  conditions  of  best 
agreement  of  calculated  pressures  Pk(x)  and  pressures  measured  along  the 
electrode  wall,  and  from  the  correspondence  of  the  value  of  uxBxhk  to  the 
electromotive  force  measured  on  a  single  remote  electrode  pair  at  the  channel  end. 

The  set  of  equations  is  solved  by  applying  the  program  for  strict  differential 
equations  similar  to  that  described  in  [44],  The  program  is  compiled  so  that  a 
comprehensive  calculation  could  be  performed  beginning  with  the  primary  nozzle 
throat  and  ending  with  the  secondary  nozzle  exit.  In  the  regions  beyond  the 
MHD-channel  the  equation 


A[(p+pu*)f]=p; 


is  used  instead  of  (8). 

The  calculations  are  conducted  untill  the  calculated  value  of  P«)  coincides  with  the 
measured  one  assuming  that  the  chemical  processes  are  frozen  except  for 
vibrational  relaxation.  Figs.  10,  12,  13  show  typical  calculated  gas  Parameter 
distributions  along  the  gas  dynamic  circuit.  The  vibrational  temperatures  in  the 
MHD-channel  proved  to  be  close  to  the  translational  ones  due  to  an  increased  rate 
of  vibration  deactivation  in  the  presence  of  O  and  N  atoms.  The  role  of  alka  i 
atoms  in  this  process  under  the  conditions  in  the  channel  is  lnsignifica  . 
Meanwhile, 'the  chemical  gas  composition  differs  slightly  from  that  in  equilibrium 
along  almost  the  whole  channel  length  (Fig.  12). 

The  high  temperature  and  increased  pressure  at  the  channel  inlet  are  caused  by 
shunting  the  Hall  currents  which  is  modelled  in  the  program  as  energy  supp  y 
gas  without  its  acceleration. 

The  alkali  atoms  in  the  secondary  nozzle  increase  somewhat  the  9as  t®f  per^r®  !j 
the  expense  of  reduction  in  Tvl  (Fig.  13)  which  causes  some  variations  in  the  flow 

Mach  number  as  well. 

The  extent  of  correlation  of  the  solutions  of  one-dimensional  equations  to  real  flow 
parameters  depends  on  the  level  of  energy  supply  uniformity  across  the  channel, 
on  turbulent  mixing  processes  in  the  channel  etc.  These  aspects  have  not  b 
adequately  studied  at  present,,  although  measured  pressures  P0  in  the  test  se 
indicate  the  presence  of  a  uniform  flow  core. 

The  comnarison  of  flow  velocities  and  the  chemical  gas  composition  at  the 
MHD-channel°exit  calculated  by  a  standard  technique  and  taking  into  account 
r»eq»Uibrium  in  the  channel,  shorn  that 
example,  the  flow  velocity  for  an  equilibrium  flow  in  the  MHD-channel  is  5100  m/ 
SH*  for  a  nonequilibrium  is  4895  m/s.  The  temperatures  are  5700  K  and 
5720  K,  the  densities  1.11-_210  kg/m3  and  1.14-10  kg/m  ,  respectively. 

The  degree  of  thermochemical  process  nonequilibrium  for  flows  over  models  in 
different1  wind  tunnels  was  evaluated  by  using  the  analog  of  the  Damkohler 
number  [15]  in  the  form: 
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cd;  =  div(piU) 

1  Poo^* 

If  ¥<1  the  flow  is  frozen.  The  case  of  T>1  corresponds  to  an  equilibrium  flow.  At 
T~1  it'  is  necessary  to  account  for  kinetics  processes.  In  the  present  study,  the 
estimations  are  performed  for  nitrogen  since  it  is  the  mam  test  gas  componen  . 
vFf  we  use  the  expression  from  [15]  in  the  form: 


=  Kh 


exp  - 


113200 


:  2r(lzi)  m2T; 

(Y  + 1) 


:  4.8  •  1018C)v 


^  teSt  "10d6lS  Wele  ^^y 

SSEjSS?  coppe?  surfaces 

vary  considerably  when  oxide  film  forms. 

The  models  were  Upped  .0  measure  pre.™  ZZ 

are  given  in  Figs.  15, 16. 

.  ,  ,  ,  -r r a  T  in9  model  pressure  distributions  were  measured  oy 

error  not  exceeding  tvthd- acceleration  wind  tunnel,  pressure 

dSiSSs^t  rs-x  » ■ 

computer.  For  these  measurements,  error  was  not  also  m  excess  of  3/,. 

SSaSSSSSSSK a 

to  correspond  to  th. 

shock  wave  front  position.  f 

£££ J£KT«2,.SS‘"  "fuciri 

In  single  tests  earned  out To' d?  tu?  1  special 
d^cT™  S«T»d  mSulaetured  since  it  vms  impossible  to  use  a  standard 
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beH.arteo°f«bi^n.r  dtjSLTrf  tohSTe  *S^  «<» 

IringtW  X~5890A  when  passing  tough  N.  vapor  presen,  to  the  flow. 

The  device  is  shown  schematically  in  Fig.  17. 

r  4.1,  to  7300A,  the  emission  line  width 

^,5  Pto/s! U  duration  .  -  — •  ^ 

i  +m°/  The  radiation  wavelength  was  checked  during  tests. 

5=S2bsi=s*“-  “  “ 
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numerical  techniques  used  in  the  present  study 

Qualitative  characterisations  of 

litii 

facilities  and  all  test  models  in  the  numerical  analysis. 

Mnmpriral  simulation  of  Chemically-nonequilibrium  hypersoms 
^Tfinw.  hv  using  th°  potions  for  a  thmvisrous  shock  layer 

■ircT  omiatinm;  RSI  includes  the  equations  of  continuity, 
momentum  projected  onto  coordinate  axes,  fhS 

S'Gte*^ 

—  »  -  «  - 

body  surface  (rj  =  0) .  .  f 

ippitnfiS  co^r^ 

parabolic  equations.  _  . 

A  chemically-nonequilibrium  gM  model  was 

solve  the  VSL  equations.  In  mathemahcrd  modeUmg  o  ttiemocn  me 

jSl^b^ 

Y.B.  Zeldovich: 

02  +  X  +  5.12  ew  <r>  20 +  X, 

N2  +  X  +  9.76  ew  o  2N  +  X, 

NO  +  X  +  6.49  ew  o  N  +  O  +  X, 

O  +  N2  +  3.27  ew  o  N  +  NO, 

O  +  NO  +  1. 37  ew  o  N  +  02 , 

02+N2+l-9ew  <->  2NO, 

where  X  is  a  catalytic  particle 

component  The  constants  ol  eb™'*"™  “  qlve„  theimochemical  model 

here  by  using  the relations  «“»™»d«d^  I«l- Tbeg^  iM^ousiy  to  a 

Uiertnodynamic'equbihrium  state  with  translational  degrees  ot  freedom. 

"  shouThe  noted  that  the c^ues  -  £  “  SrtS 

S»tr,cS“5  wa?£«totoP«“  ™n.  ta”ES— .  -  -  results. 
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A  library  of  universal  programs  designed  to  solve  parabolic  nonlinear  differential 
equations  in  an  undetailed  form  was  used  for  numerical  integration  of  the  VSL 
equations  describing  hypersonic  chemicallynonequilibrium  gas  flows  over  blunt 
bodies  [471.  The  two-point  second-order  Keller  scheme  was  applied  for  the  solution 
of  the  VSL  equations  to  approximate  the  derivatives  in  the  direction  normal  to  the 
body  surface  [48].  The  advantage  of  this  scheme  is  that  the  problem  of  boundary 
layer  approximation  can  be  easily  solved  irrespective  of  grid  step  ratios.  The  Keller 
scheme's  nonmonotony  was  overcome  by  applying  an  adaptive  grid. 

An  absolutely  stable  and  fully  implicit  second-order  scheme  was  used  to 
approximate  longitudinal  derivatives  in  the  VSL  equations.  This  scheme  features 
good  dispersion  properties  and  damps  short-wave  oscillations  of  grid  functions 
rather  effectively. 

The  qeneral  procedure  of  solving  two-  and  three-dimensional  problems  consisted  of 
two  stages.  The  first  one  included  the  solution  of  the  self-similar  VSL  equations  on 
the  flow  stagnation  line  by  using  an  adaptive  grid  for  the  variable  t|.  A  nonuniform 
qrid  constructed  in  the  solution  of  the  self-similar  problem  was  fixed  in  most 
calculations  to  use  it  at  the  second  stage  for  marching  solution  of  parabolic  two- 
and  three-dimensional  VSL  equations. 

The  Newton-Raphson  method  [49]  was  used  for  numerical  solution  of  nonlinear 
qrid  equations.  The  correction  vector  for  each  iteration  of  this  method  was 
identified  by  solving  linearized  equations.  The  vector  factonng  was  applied  o 
solve  linear  equations  when  dealing  with  one-dimensional  (stagnation  line)  and 
two-dimensional  (separation  line)  VSL  equations. 

As  for  the  boundary  conditions  on  the  body  surface,  the  wall  temperature  in  all 
cases  was  assumed  to  be  Tw=300  K  while  the  wall  itself  was  assumed  to  be 
absolutely  noncatalytic. 

Mnmpriral  simulation  of  hypersonic  gas  flows  over  blunt  bodies 
using  full  navjerzstokes  equations  and  viscous  shock  layer  equations 

Perfect  supersonic  perfect  gas  flows  over  a  sphere,  a  cone  and  a  wedge  featuring 
heat-insulated  surfaces  were  calculated  by  using  the  shock-captunng  method  based 
on  full  Navier-Stokes  equations  [50],  [51]. 

The  Navier-Stokes  equations  in  a  strictly  conservative  form  were  integrated 
numerically  over  a  finite  region  limited  by  external  boundary,  body  surface  plan 
of  symmetry,  and  exit  boundary.  The  impermeability  and  nonslip  conditions,  as 
well™  the  condition  of  heat-insulated  body  surface  were  posed  on  ^he  internal 
boundary.  The  dependence  of  the  pressure  gradient  normal  to  the  surface  on  the 
second  Velocity  derivative  obtained  asymptotically  from 

was  used  as  an  additional  boundary  condition  on  the  body  surface.  This  was 
required  to  represent  the  grid  equations  in  full.  The  ccnditions  of  evennes 
oddness  for  sought  functions  were  posed  along  the  line  of  symmetry, 
conditions  of  undisturbed  flows  were  used  on  the  external  boundary  and  the 
conditions  of  extrapolation  of  dependent  variables  on  the  outflow  boundary. 

The  finite-difference  scheme  for  equations  written  in  the  form  of  conservation  laws 
was  based  on  an  integro-interpolation  method  ensunng  the  fulfilment  of 
conservation  laws  (for  the  Euler  equations  even  in  the  presence  of  discontmuou 
solutions).  The  convective  flow  vector  component  at  half -“teg^°£  Se 
armroximated  by  using  a  monotonic  second-order  scheme  [52],  [53]  ihe  k 
method  for  an  approximate  solution  of  the  problem  of  the  break-down  of  arbitrary 
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noder'lS^I^Th^dtffusive1  Vtoear 

equations  (including  mixed  derivatives)  included  13  points. 

a  moHified  Newton  method  was  used  to  solve  a  nonlinear  set  of  grid  equations^ 
Th“  Jacobi  matrix  was  formed  in  each  iteration  by  applying  a  procedure  of  finite 
teKmenS  S  residual  vector  in  terms  of  the  vector  of  sought  gnd  funchons. 

The  linear  algebraic  equations  resulted  from  each  nonlinear  «erationr vere  solv, ad 

totSrg!  S  “dues  were  renumbered  by  using  the  generated  nested 

dissection  method  [55]. 

Reynolds  number. 
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TEST  DATA  AND  THEIR  ANALYSIS 


The  gasdynamic  flow  parameters 
in  the  following  tables: 


used  to  investigate  model  flows  are  summarized 

Table  1 


Wind  tunnel 

n 

Q 

EM 

Lu 

VAT-102 

Nitrogen 

lip 

7.63 

Arc-heating  wind 
tunnel 

7.97 

MHD -acceleration 
wind  tunnel 

1^^ 

B 

Table  2 


Gas  composition  (mole  fractions) 

- - - - - - 

Wind  tunnel 

N2 

o2 

NO 

N 

o 

K,Na 

VAT- 102 

1.0 

- 

- 

- 

- 

Arc -heating 

wind  tunnel 

7.4510"1 

1.66-10“ 1 

4.2- 10“2 

1.6-10“8 

3.710"2 

10“2 

0 

MHD- acceleration 

5.88- 10“ 1 

3.09- 10- 4 

5.04- 10"3 

6.9- 10“2 

3.29- 10“ 1 

10“2 

"5  2WK 

the  main  tests  were  conducted  with  nickel-coated  models.  As  a  result, 
surfaces  remained  clean  throughout  the  balance  of  testing. 

In  tests,  model  surface  pressure  distributions  were  obtained,  as  well  as  the 
information  about  shock  wave  shapes  and  stand-off  distances. 

Fias  18  -  20  show  surface  pressure  distributions  for  a  hemisphere, a  cone,  and  a 
wedqe  Pressure  values  are  related  to  respective  dynamic  pressures  pu  for  toe 
used.  For 

d“dTu^el1bvr3th^Se CrSsV pressure  d^rtbuUo^ 

"  Sr 

corresponding  to  the  solution  of  the  Navier-Stokes  equation  for  a  perfect  gas. 
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The  flow  visualization  data  are  presented  on  photos  and  video  frames  showing 
natural  gas  glow  near  models,  on  shadow  pictures  for  flow  fields,  as  well  as  on  all 
these  pictures  after  computer  processing. 

A  photo  of  gas  glow  over  a  model  obtained  in  wind  tunnel  VAT- 102  by  using  glow 
discharge  and  its  respective  video  frame  are  shown  in  Fig.  21,  and  photos  and 
video  frames  of  natural  gas  glow  over  models  in  the  hypersonic  arc-heating  and 
MHD-acceleration  wind  tunnels  in  Figs.  22  to  25. 


Fig.  26  displays  typical  computer-processed  video  pictures  for  gas  glow  over 
models  along  nitrogen  molecule  lines. 

For  ease  of  comparison  with  theoretical  results,  the  visualization  data  for  a 
hemisphere  were  reduced  to  dependences  of  relative  shock  wave  stand-off 
distances  normal  to  the  model  surface  on  relative  distances  along  its  surface  from 
the  critical  point.  Cone  and  wedge  data  were  reduced  to  dependences  of  shock 
stand-off  distances  normal  to  their  surfaces  on  distances  along  the  axis  from  the 
model  nose.  Some  examples  of  this  data  reduction  are  given  in  Figs.  27  to  29.  Later 
on  the  data  of  only  those  runs  were  used  when  the  influence  of  model  installation 
inaccuracy  was  excluded,  i.e.,  when  shock  stand-off  distances  for  the  upper  and  the 
lower  surfaces  coincided. 

Fias  30  to  33  qive  the  summarized  flow  visualization  data  for  all  models  tested  and 
the  calculated  results  obtained  using  VSL  theory  and  the  Navier-Stokes  equations 
based  on  experimental  flow  parameters  for  each  particular  run. 

Comparison  of  experimental  and  theoretical  pressure  distributions  over  a 
hemisphere  (Fig.  18)  shows  that  both  for  a  perfect  gas  and  for  a  gas  with  Physical 
and  chemical  processes  the  experimental  data  are  in  better  agreement  wi 
distribution  calculated  using  the  Navier-Stokes  equations  for  a  perfect  gas  than 
with  that  obtained  by  VSL  theory.  The  discrepancy  becomes  appreciable  at 
considerable  distances  from  the  critical  point  which  can  be  attributed  to  the  fact 
that  the  Navier-Stokes  equations  describe  flows  over  a  hemisphere  m  full  There 
are  no  terms  in  the  VSL  equations  which  are  responsible  for  convective  transfer 
and  diffusion,  and  the  surface  pressure  distribution  is  assumed  to  he?OTI^°^e 
to  a  modified  Newton  distribution.  The  VSL  theory  is  valid  only  in  the  vicinity  of 
the  forward  critical  point.  This  figure  also  shows  that  physico-chemical  processes  in 
gas  have  slight  influence  on  the  hemisphere  pressure  distribution.  As  for  a  cone 
and  a  wedge,  the  discrepancy  for  them  is  more  significant  (Figs.  19,  2^.  For  a 
cone  initial  pressure  values  correspond  to  theoretical  ones  for  a  Perfect  9a 
(curve  1).  TtJpressure  along  the  cone  generatrix  decreases  quicker  than  it' follows 
from  the  Navier-Stokes  equations.  Respective  pressure  distribution  calculated  by 
VSL  theory  lies  lower  than  the  experimental  data.  The  <>“arepancy 
theoretical  and  experimental  data  is  likely  to  be  caused  by  e 
calculation  of  flows  over  a  cone  and  a  wedge  by  using  the  Navier-Stokes  equations 
?  assuming  infinite  extension  and  without  accounhng  lot  chern^l 

reactions,  and  the  calculation  by  applying  VSL  theory  proceeding  from  the  above 
considerations. 

The  above  stated  assumptions  should  be  verified  in.  the  future  by  testing  models 
with  great  relative  generatrix  lengths  and  by  comparing  the  data  obtained  wth  the 
calculations  using  the  Navier-Stokes  equations  including  physical  and  chemical 
processes. 

The  analvsis  of  the  data  on  shock  wave  shapes  and  stand-off  distances  for  a 
hemisphere  (Fig.  30)  shows  that  for  a  perfect  gas  the  experimental  points  lie  stnc :  y 
on  the  curve  Calculated  by  the  Navier-Stokes  equations.  At  the  same  time,  the 
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shock  front  positions  determined  according  to  gas  glow  boundaries  in  flow 
conditions  2  and  3  differ  greatly  both  from  each  other  and  from  the  theore 

data  The  test  points  for  flow  condition  2  are  lower  than  the  points  caicuiated  by 
th  N  vier- Stokes  equations  and  somewhat  higher  than  the  points  obtained  y 

corresponds  to  a  spherical  surface  up  to  angles  of  -50°.  After  that,  the  test  points 
also  lie  on  the  Navier-Stokes  curve. 

The  shock  wave  position  determined  from  shadow  pictures  differs  significantly 
while  at  areat  angles,  the  test  points  correspond  to  the  NS  curve,  it 

Ss  “  ££  SMLS  « 

especially  for  flow  condition  3. 

wmm mMmm 

the  nose.  . 

The  above  considerations  based^n  different 

and  so  on.  .  ,  .  , 

A  number  of  effect  were  revealed  during 

critical  point  for  a  hemisphere  in  flow  condition  3,  local  flow 

nose  etc.  It  would  be  possi  effects  the  test  models  were  calculated  within  the 
characteristics.  To  clear  up  these  eff  ,  vibrational  temperatures  behind 

framework  of  the  VSL  “^T.empSS  dSes  and  particle 

sissL”  —  — *  ““VFS'stoS1  ar“ 'low 

conditions.  The  results  of  these  calculations  are  given  m  Figs.  34  to  4  . 

It  is  seen  that  there  are  drastic  temperature  and  density  gradient, »  near  the 

shS^e^Soff  dXcflcSrdi^  to  t“cuTai°  tie  temperature 
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factor  TwX  sb  0.03  has  an  influence  only  at  a  very  small  distance.  Similar  effects 

/To 

also  take  place  for  a  cone  and  for  a  wedge. 

Figs.  36  to  38  show  concentration  variations  of  monatomic  and  molecular 
components  in  a  shock  layer  over  a  hemisphere  at  different  distances  from  its 
surface.  It  is  seen  that  an  additional  dissociation  of  nitrogen  and  oxygen  takes 
place  in  a  shock  layer,  while  recombination  processes  begin  to  take  effect  only  in 
the  close  vicinity  of  a  cold  wall.  Of  importance  is  a  drastic  variation  in  NO 
concentration  along  streamlines.  Similar  results  are  obtained  in  shock  tunnel  tests. 
As  stated  above,  the  excitation  of  electron  states  of  NO  molecules  causes  an 
intensive  ultraviolet  glow  from  this  flow  region  [6], 

On  the  other  hand,  monatomic  components  of  N  and  O  would  result  in  a  certain 
influence  of  the  catalycity  level  on  shock  stand-off  distances  [16].  In  the  given 
investigations,  the  influence  of  catalycity  level  on  shock  stand-off  distances  is  not 
considered.  Moreover,  nickel  featuring  weak  catalytic  properties  is  chosen 
deliberately  for  model  surfaces. 

A  possible  difference  of  the  glow  front  position  from  the  shock  front  may  be 
explained  by  two  mechanisms:  by  removal  of  glow  from  the  shock  layer  region 
basically  in  the  form  of  resonance  radiation  and  by  emission  of  high-energy 
electrons  from  a  high-temperature  region  immediately  adjacent  to  the  shock  wave 
front  These  two  phenomena  can  result  in  an  excitation  of  electron  states  o 
impinging  molecules  and  atoms  of  air  and  seeds  which  may  be  attributed  to  the 
presence,  in  flow  condition  3,  of  a  sharp  front  boundary  of  the  glow  region  at  the 
critical  point  of  a  hemisphere  and  its  spherical  character  up  to  anges 
corresponding  to  the  transition  of  a  subsonic  region  to  a  supersonic  one.  The 
effects  of  radiation  removal  are  characteristic  of  strong  shock  waves  caused,  for 
example,  by  nuclear  explosion  [56]  or  of  meteors  entering  the  atmosphere  [57], 

The  situation  close  to  the  flow  condition  considered  here  is  discussed  in  [58], 
Based  on  numerical  solution  of  the  Biberman  equations  for  radiation  diffusion  in 
resonance  lines  it  is  shown  that  the  effect  observed  may  be  explained  by  radiation 
transfer  in  resonance  lines  Cu  and  Na  with  their  subsequent  reradiation  to  the 
region  ahead  of  the  shock  wave.  In  flow  condition  3,  the  above  effects  take  place 
not  only  for  seeds  having  a  small  excitation  potential  but  also  for  nitrogen  and 
oxygen.  The  analysis  of  these  processes  is  beyond  the  scope  of  this  study. 

The  numerical  estimates  of  characteristic  electron  free-path  lengths  as  applied  to 
the  conditions  of  emission  from  a  high-temperature  gas  show  that  they  are  ot 
~2  mm,  i.e.,  they  are  close  to  experimental  glow  zone  dimensions.  It  is  known,  on 
the  other  hand,  that  a  body  moving  with  velocities  of  ~5+6  km/s  can  gam  electric 
charge  of  a  order  of  several  V  [59]  which  also  contributes  to  electron  emission. 

In  conclusion,  consider  the  following  point.  An  obligatory  condition  of  the 
MHD-acceleration  wind  tunnel  operation  is  the  presence  of  alkali  seeds  in  the 
flow.  In  the  given  tests,  seeds  were  used  in  single  runs  and  m  the  arc-heating  wind 
tunnel.  Although  their  small 

fraction  (~1%  of  weight)  could  not,  naturally,  change  thermodynamic  gas 
parameters  (pressure,  enthalpy,  heat  conductivity  etc.)  they  could  influence  rates  of 
relaxation  processes.  In  order  to  clarify  this  point,  methods  of  mathematical 
modelling  were  applied  to  investigate  translational  and  vibrational  temperature 
distributions,  as  well  as  air  component  concentrations  behind  a  normal  shock  wave 
and  shock  waves  being  at  angles  of  51°  and  38°  to  the  flow.  The  last  condition 
corresponded  to  test  shock  wave  positions  for  a  cone  and  for  a  wedge,  respective  y. 


23 


CONTRACT  CPS-93-4044 


The  calculation  was  performed  for  the  same  free-stream  flow  parameters  without 
Na  and  K  and  with  them.  The  concentration  was  assumed  to  be  equal  to  -1%  with 
respect  to  gas  mass  flow  rate.  The  calculated  results  are  given  in  Figs.  43F  44  . 

It  is  seen  that  in  the  case  of  a  normal  shock  wave,  the  seeds  exert  a  very  small 
influence  on  thermodynamic  processes.  The  attainment  of  thermodynamic 
equilibrium  at  given  parameters  required  ~10~4  s.  The  influence  of  seeds  on 
relaxation  of  vibrational  degrees  of  freedom  is  most  appreciable  behind  shock 
waves  near  a  wedge  and  a  cone.  But  in  this  case,  their  influence  on 
thermodynamic  processes  (P,  T)  is  insignificant  (Fig.  44).  Some  dispacement  of  the 
nitrogen  vibrational  temperature  maximum  relative  to  translational  one  can  result 
in  displacement  of  a  radiation  front  relative  to  a  shock  wave  since  electron 
temperature  is  governed  here  by  vibrational  temperature.  This  phenomenon  can  be 
a  cause  of  radiation  delay  in  flow  condition  3  near  a  cone  nose  (Fig.  32)  and  it  is 
likely  to  be  attributed  to  the  coincidence  of  the  shock  wave  determined  by  the 
shadow  technique  with  that  calculated  using  the  Navier- Stokes  equations  for  a 
perfect  gas.  Indeed,  the  shock  wave  position  for  a  fully  frozen  gas  would 
correspond  to  a  flow  with  y=1.48.  Due  to  relaxation  of  vibrational  degrees  of 
freedom  the  quantity  y  decreases  and  approaches  a  value  close  to  1.40.  The  fact 
that  the  shock  wave  stand-off  distance  in  the  vicinity  of  cone  nose  is  greater  than 
that  followed  from  the  Navier-Stokes  equations  confirms  this  phenomenon  since  in 
this  case  the  flow  in  the  cone  nose  vicinity  is  close  to  a  frozen  one. 
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CONCLUSION 

Flows  over  simple  bodies  (a  hemisphere,  a  cone,  a  wedge)  have  been  investigated 
in  wind  tunnels  of  different  types  at  close  similarity  parameters  for  Mach  and 
Reynolds  numbers  but  at  essentially  different  flow  velocities.  The  tests  were  earned 
out  in  the  vacuum  wind  tunnel  VAT- 102  for  the  conditions  of  a  perfect  gas 
(u  =  790  m/s,  M  =  7.63,  Re0  =  240,  H0  =0.32  MJ/kg),  in  the  arc-heating  wind  tunne 
in  a  gas  flow  featuring  excited  vibrational  degrees  of  freedom  and  a  sma 
dissociation  level  (u  =  2700  m/s,  M=7.97,  Re0=150,  H0  =  5.2  MJ/kg)  and  m  the 
MHD-air  flow  acceleration  where  the  influence  of  physical  and  chemical  processes 
can  be  significant  (u  =  5950  m/s,  M=7.23,  Re0=  135,  H0=29MJ/kg). 

The  matched  operation  conditions  of  the  facilities  used  and  the  means  of  their 
checking  were  optimized  during  the  tests  by  applying  the  technique  of 
determining  both  gasdynamic  (M,  Re,  Re0,  T,  P,  U,  Ho,  H)  and  physical  and 
chemical  (y,  T0,  TVf,  Ttr,  Q)  free-stream  flow  parameters.  Particular  attention  was 
given  to  the  influence  of  nonequilibrium  processes  in  a  gas  accelerated  in  the 
MHD-channel  of  the  hypervelocity  wind  tunnel  and  in  the  secondary 
shown  that  nonequilibrium  effects  are  insignificant  in  the  MHD-channel.  The  flow 
in  the  secondary  nozzle  is  closer  to  a  frozen  one;  alkali  seeds,  being  in  air  flow, 
approach  the  thermodynamic  gas  to  an  equilibrium  state  mainly  y  acce  era  i  g 
deactivation  of  vibrational  degrees  of  freedom  of  nitrogen. 

Different  means  of  flow  visualization  are  devised  and  tried.  They  are  based  on 
recording  of  flow  glow  near  models  excited  by  glow  discharge  electrons  (  ) 

or  natural  flow  glow  (arc-heating  and  MHD-acceleration  wind  tunnels).  Both 
neutral  and  interference  light  filters  were  used  in  the  MHD-acceleration  wind 
tunnel  A  shadow  device  has  been  developed  based  on  the  effect  of  abnon™* 
dispersion  of  Na  atom  resonant  radiation  with  the  application  of  an  impulse  las 
of  variable  radiation  frequency  as  a  light  source.  This  made  it  possible  o  m 
1  SnsSy  by  ,  (actor  of  10  and  to  exclude  distortions  due  to  shoclt 

wave  front  glow.  : 

Pressure  distributions  over  model  surfaces  were  measured  in  all  wind  tunnels  used 
and  compared  with  calculations  performed  for  test  conditions  by  applying 
theory  and  the  Navier-Stokes  equations . 

As  for  a  hemisphere,  the  influence  of  the  flow  velocity  on  the  pressure  distribution 
its  surfaCe  has  not  revealed  and  in  all  cases  it  was  in  good  agreement  wi 
Zl  c,lc„“£d  tSg Navier-Stokes  ted  by  the  VSL  theory  and  Navier-Stokes 

equations. 

Thus  a  number  of  new  effects  have  been  revealed  during  the  investigation  carried 
St  which  cannot  be  explained  in  the  framework  of  integral  characteristics,  tt 
seems  to  be  possible  that  they  depend  on  local  shock  layer  parameters.  To  verity 
thfs  suacrestion  they  were  calculated  by  applying  the  VSL  theory;  the  influence  of 
alkali  seeds  on 'relaxation  rate  was  analyzed.  Some  fields  of  gasdynanuc 

sets  Mo 

addition^the  available  high-temperature  region  is  a  source  of  rapid  electrons.  Both 
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these  factors  can  cause  radiation  removal  for  a  distance  behind  the  shock  wave 
front. 

It  is  shown  that  seeds  of  -1%  behind  a  normal  shock  wave  do  not  influence  the 
chemical  reaction  rate  and  vibrational  temperature  relaxation,  and  behind  an 
oblique  shock  wave  (cone,  wedge)  its  influence  is  very  weak  (~2  to  3%).  But  behind 
an  oblique  shock  wave  the  vibrational  temperature  maximum  is  father  away  from 
the  shock  wave  front  because  of  a  greater  gas  velocity  which  can  shift  the  radiation 
intensity  peak  closer  to  the  body  surface. 

The  results  obtained  also  show  that  it  is  important  to  choose  a  required 
gasdynamic  model.  It  is  follows  from  the  comparison  of  theoretical  and 
experimental  data  that  the  flow  model  based  on  the  Navier-Stokes  equations  is  in 
better  agreement  with  the  test  results  than  that  calculated  using  the  VSL  theory. 
Therefore,  the  multitemperature  kinetics  of  chemical  reactions  in  air  should  be 
considered  within  the  framework  of  the  Navier-Stokes  equation  solution.  The 
reproduction  of  a  stationary  flow  field  over  a  model  and  that  repeated  from  run  to 
run  allows  the  statement  of  a  problem  concerning  the  evaluation  of  a  detailed 
structure  of  the  shock  layer,  i.e.,  chemical  air  component  distribution  in  the  layer, 
both  equilibrium  and  nonequilibrium  radiation  intensity  and  spectrum. 

It  seems  to  be  reasonable  to  evaluate  the  influence  of  the  model  surface  catalycity 
level  on  the  shock  wave  stand-off  distance  and  the  flow  character.  The  capabilities 
of  available  hypervelocity  wind  tunnels  and  other  hypersonic  wind  tunnels  of 
TsAGI  make  it  possible  to  carry  on  similar  experimental  investigations  of  greater 
models  at  considerably  higher  Mach  numbers  M— 15  to  20  but  at  slightly  less 
Reynolds  numbers  Res  =  50  to  100. 
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accelerator  channel  width 

magnetic  field  induction 

channel  drag  coefficient 

concentration  of  i-p  articles  (mole  fractions) 

diameter 

hydraulic  diameter 
energy 

vibrational  energy 
vibrational  energy  per  i-particle 
electric  field  strength 

energy  required  to  generate  i-component 

channel  cross-section  area 

current  density 

channel  height 

enthalpy 

current 

function  of  generation  of  i-particles 

characteristic  size,  channel  length 

molecular  weight 

gas  mass  flow  rate 

Mach  number 

pressure 

heat  flux  related  to  unit  volume 

coefficient  responsible  for  electrode  current  nonuniformity  [30] 
radius 

universal  gas  constant 
Reynolds  number 

Reynolds  number  based  on  stagnation  parameters 

channel  (electrode)  sectionning  step 

temperature 

velocity 

measured  electrode-to-electrode  voltage 
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near-electrode  voltage  drop 
coordinate  along  flow 

coordinate  across  flow  normally  to  applied  magnetic  field  B 

energy  transfer  rate  to  vibrations  of  i-particles 

total  energy  fraction  supplied  to  gas 

total  energy  fraction  tranferred  to  N2  molecules 

Hall  parameter 

specific  heat  ratio 

shock  wave  stand-off  distance 

viscosity 

density 

electric  conductivity 
characteristic  time 

Damkeller  number 
circular  frequency 

source  term  in  chemical  reaction  for  i-particles 


Subsripts 
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chem  chemical 
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diss  dissociation 
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e  electrode 

eq  equilibrium 

fr  frozen 

g  gas 

H  Hall 

i=  1,  2,  3,  4,  5,  6  reffer  to  N2,  02,  NO,  O,  N,  Na,  respectively 

k  electrode  number 

m  model 

n  full-scale 

o  total  (stagnation) 

Tr  translational 

Y  vibrational 

vtr  vibrational-translational 

VV'  vibrational-vibrational 
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Fig. 2.  Pressure  field  P'o  VAT  —  102. 
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Fig.,.  Schematic  of  hypecveioci.y  MHD-  ««■«““  wind  tunnei 
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\  —  arc  —  heater; 

2  —  mixing  chamber; 

3  —  metering  device; 

4  -  primary  supersonic  nozzle, 

5  —  MHD  —  accelerator; 

0  _  secondary  nozzle, 

7  —  test  section 


4 


I 


Pig. 8.  Hypervelocity  MHD-gas 


Fig.  9.  Hypervelocity  MHD  —  gas  acceleration  wind  tunnel.  General  view. 
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Fig.  10.  Electric  and  gasdynamic  parameter  distributions  along  the 
accelerator  channel 
1—  temperature;  T  •  10  ,K 

2  -  velocity;  u-10 

3  -  static  pressure;  P-10'5,Pa 

4  —  density;  P  ^  2  ^ 

5 -electrode  current;  M0  ,A 

fi  —  electrode  voltage;  V-10"3,1^ 


Fig. 11.  Pressure  field  P'o 

without  MHD- acceleration;  2  -  with  MHD  -  acceleration 
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Fig.  14.  Model  drawings 


Fig.  17.  Schematic  of  a  modernized  shadow  device 

1  —  Chramatic  CMX  —  4  laser; 

2  —  light  guide; 

3r8  —  shadow  device; 

4  —  rotatiall  mirror; 

5  —  test  section; 

6  —  test  model; 

7  —  adjustable  rotating  mirror; 

9  —  Fucoult  knife; 

10  —  light  filters; 

1 1  r  1 2  —  scintific  — grade  video  camera  and  video 


Fig.  18.  Pressure  distributions  over  a  hemisphere 
VAT— 102;  2  —  arc  — heater;  3  —  MHD  —  accelerator 
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Fig.  21.  Hypersonic  wind  tunnel  VAT  —  102. 

a)  photo 

b)  video  frame. 


Fig.  22.  Arc-heater.  Flow  pictures. 


Fig.  24.  MHD- accelerator  (pictures  on  the  left  side  are  obtained  by  using  an 
interference  light  filter  and  those  on  the  right  side  by  a  neutral  light 
filter) 

a)  photo 

b)  video  frame 
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Fig. 29.  Flow  glow  zone  stand-off  distances  (a)  and  shock  wave 
stand-offs  (shadow  pictures)  for  a  cone. 


Fig. 30.  Flow  glow  zone  stand-off  distances  for  a  hemisphere  (test  data 
are  shown  by  points  and  calculated  data  by  lines) 

1  —  VAT— 102;  2  —  arc  —  heater  wind  tunnel; 

3  -  MHD  — acceleration  wind  tunnel 


Fig.  31.  Shock  wave  stand-off  distances  for  a  hemisphere  obtained  from 
shadow  pictures  (test  data  are  shown  by  points  and  calculated 
data  by  lines) 

1  —  VAT— 102;  2  —  arc —heater  wind  tunnel; 

3  —  MHD  —  acceleration  wind  tunnel 


Fig. 33.  Flow  glow  zone  stand-off  distances  for  a  wedge  (test  data  are 
denoted  by  points  and  calculated  data  by  lines) 

1  —  VAT  —  102;  3  —  MHD  —  accelerator 
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ig.35.  Density  field  in  a  shock  layer  near  a  hemisphere 
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Fig.43.Variations  of  temperatures  and  concentrations  of  air 
components  and  air  with  Na  and  K  seeds  with  time  behind  a 
normal  shock  wave. 


